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CONVECTION  IN  A CLOSED  CAVITY  HEATED  FFCM  TEE  SIDE  WITH  THE 
DEPENDENCE  OF  VISCOSITY  ON  TEMPEF/TIJFE 

Ye.  L.  Tarunin  and  V.  I.  Chernatynski y 


During  the  theoretical  investigation  ot  free  thermal  convection, 
it  is  usually  assumed  that  the  parameters  which  characterize  th^ 
physical  properties  of  the  fluid  are  constants.  Actually,  however, 
these  values  depend  on  temperature  and  fresstre.  In  many  cases,  th« 
dependence  of  the  viscosity  coefficient  on  temperature  is  especially 
significant.  H igh- v iscosity  fluids  - glycerin,  cils,  certain 
petroleum  products,  and  silicone  fluids  - exhibit  a strong  dependence 
cf  viscosity  on  temperature.  Even  for  water,  in  the  temperature  range 
frem  10°C  to  100°C  the  Kinematic  viscosity  coefficient  undergoes  more 
than  a four- fold  change. 
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The  investigation  cf  the  dependence  of  viscosity  on  temperature 
entails  considerable  mathematical  difficulties;  therefore, 
theoretical  investigations  cf  convection  with  consideration  of  ‘his 
dependence  are  rare  [1-6].  The  precise  soluticr  of  the  problem  of 
stable  convection  in  an  infinite  vertical  flat  layer  was  obtained  in 
in*  The  stability  of  flov«s  originating  in  a horizontal  flat  layer 
heated  from  below  was  considered  in  [2,  3],  Eepcrt  [4]  considered 
movement  in  a convective  cell  with  free  boundaries  as  a model  of 
movement  in  the  Earth's  upper  mantle.  The  problem  of  convection  in  a 
vertical  cylindrical  container  was  solved  ir.  [5],  The  methol  of  nets 
was  used  t.  c consider  the  problem  of  convection  in  a perfect  gas  in 
[6]. 

This  report,  gives  the  results  of  the  numerical  study  of 
convection  in  horizontal  cylinders  with  square  and  round  cross 
sections  heated  from  the  side. 

1.  Formulation  of  Problem.  Main  Equations. 

When  viscosity  varies,  in  vector  form  the  equations  o*  ‘"re* 
ccnvecticn  are 
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/in  — ♦ — * I -* 

-37  + (it)*  = - — r/>  + + 2(^r)v-r 

+ (rvx  rot  u)  + (1.1) 

tr  (i.2> 

div  u *0.  (1.3) 

Here  v - velocity,  T - temperat are , p - pressure,  p0  - the  s°an 
density  of  the  fluid,  g - the  acceleration  of  the  force  of  grav’*yf 
and  x - the  coefficients  of  thermal  expansion  arl  heat  conductivity, 
and  k - the  unit  vector,  directed  upward. 

Coefficients  p and  x are  considered  to  he  constant  in  equations 
(1.1)- (1.2),  while  the  kinenatic  viscosity  coefficient  v is  assumed 
to  be  knowr  from  the  temperature  function.  Ihe  linear  dependence  of 
viscosity  on  temperature  is  considered  in  this  study: 

v = v'o(I — uT).  (1.4) 

where  v0  is  the  mean  viscosity,  which  corresponds  to  the  n®in 
teirperature  taken  as  the  computing  origin. 

Two  problems  are  considered  in  this  study:  a)  convection  in  an 
infinite  horizontal  cylinder  with  a square  cress  section;  H) 
convection  in  a cylinder  with  a round  cross  section. 


We  will  formulate  the  {rc'clem  of  plane  convective  movement  ecr 


■57  + (vv) v = ---  \p  + Av  + 2(vvr)'t'  + 

+ (vvX  rot  v)  4-  g?Tk,  (1.1) 

%-  -\-v?T=  />r,  (1.2) 

cl i \ t»  = 0 . (1.3) 

Here  v - velocity,  T - temperature,  p - pressure,  p0  - the  near, 
density  of  the  fluid,  g - the  acceleration  of  the  force  of  gravity, 
and  x - the  coefficient  cf  thermal  expansion  ard  heat  conductivity, 
and  k - the  unit  vector,  directed  upward. 

Coefficients  (3  and  x are  considered  to  he  constant  in  equations 
(1-1)  — (1.2),  while  the  kinematic  viscosity  coefficient  v is  assumed 
to  be  knowr  from  the  temperature  function.  Ihe  linear  dependence  of 
viscosity  on  temperature  is  considered  in  this  study: 

v = v0(l—  uT),  (1.4) 

where  vn  is  the  mean  viscosity,  which  corresponds  to  the  m®an 
temperature  taken  as  the  computing  origin. 

Two  problems  are  considered  in  this  study:  a)  convection  in  an 
infinite  horizontal  cylinder  with  a square  cress  section;  h) 
convection  in  a cylinder  with  a round  cross  section. 


He  will  formulate  the  problem  of  plane  convective  movement  for 
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cavity  with  a square  cress  section  with  side  d (Fig.  1).  All  of  the 
boundaries  of  the  region  are  assumed  to  be  solid.  Constant 
temperatures  of  ±0/2  are  naintained  cn  the  vertical  boundaries,  while 
the  horizontal  boundaries  are  assuired  tc  be  thermally  insulated. 


f 


Earameter  7 is  related  to  tie  value  of  the  viscosity  drop  n by  the 

ratio 


P 


'max 

vmin 


2 + 7 

2 — 7 


(1.9) 


The  boundary  conditions  are: 


at  x = 0:  + = ^ = 0,  T^-~- 

at  *-l;  •>=  ^-==0,  r=  + 4-l  CV.K’) 

. r\  ^ 7*  n 

at  y = 0,  y = 1:  V = 57  = -57  ~ °’ 


/■ 


This  problem  was  solved  by  the  finite  differences  method;  an 
explicit  system  with  the  substitution  of  the  central  differences  fo 
the  spatial  derivatives  was  used.  Tte  procedure  of  the  solution  was 
described  earlier  [1,  9].  Tie  time  interval  varied,  calculated  from 
the  formula : 


Af  - 


*1 


(8  + IO 


(1.11) 


where  T is  the  maximum  absolute  value  of  the  current  functions  a* 
the  nodes  of  the  network. 
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The  computat i oris  were  conducted  on  an  Aragats  computer  or.  16x16 
and  21x21  networks.  The  valves  cf  tie  Piandtl  Timber  were  p = ^ an-1  P 
= 5.  The  viscosity  heterogeneity  parameters  selected  were  y = C,  0.6, 
C.75,  1.0,  1.25  ar.d  1.627  (the  last  value  corresponds  to  a ten-fold 

change  in  viscosity:  p = 10).  The  maximum  value  cf  the  Grashof  number 
is  G - 50» 1 0 3 . 


Now  we  will  consider  the  fcrmulaticn  cf  the  problem  for  a round 
horizontal  cylinder  wj<-h  radius  F.  The  tourcaries  are  heated 
according  to  the  cosine  law,  sc  that  tie  maximum  temperature  ft 
corresponds  to  angle  $ = C°,  and  *he  minimus  6 - to  the  value  «'  - 
180°.  In  cylindrical  coordinates  and  dimensionless  form,  the 
equations  of  plane  convective  movement  are: 
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or 

1 / fly 

ar 

&b  ,)T\ 

ft  ' • 

(1.14) 
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Here  the  Prandtl  and  Grashof  numbers  were  determined  crom  the 
mean  viscosity  v0i 
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The  current  function  was  introduced  in  he  usual  manner: 


V*~!F- 


(1.16) 


F , F 2 /d  0 , vo  and  9 are  the  units  of  distance,  time,  current  function, 
and  temperature.  It  should  le  pointed  out  that  the  total  temperature 
difference  in  this  problem  is  29,  which  causes  the  sense  o r para  me4- ->r 
7 to  change  somewhat;  tie  value  7 = 0.818  new  corresponds  to  a 
ten-fcld  drop  in  viscosity. 


The  boundary  conditions  of  the  problem  a: 


at  r 1 : 'j*  = 0,  4-r  = 0,  T — cos  *>. 


(1.17) 


The  method  of  nets  was  used  tc  solve  problems  ( 1 . 1 2) - ( 1 . 1 7)  . ? 
network  which  was  uniform  with  respect  to  tie  radius  and  angle  with 
ten  steps  on  the  radius  and  16  steps  on  the  angle  was  used  in  tbo 
(tain  calculations.  A network  which  was  irregular  over  the  radius 
(thickening  at  the  edge)  was  used  in  seme  ci  the  calculations.  ”hn 
order  of  approximation  cf  tie  equations  and  bcurdary  conditions  was 
equal  to  0 (h-2)  , where  h is  *-he  maximum  distance  between  the  nod-s  of 
the  network.  Symmetrical  differences  were  used  to  replace  the  spatial 
derivatives.  The  indeterminanoy  of  certain  expressions  in  the 
differential  equations  at  r — >0  was  revealed  by  the  l*H$oital  rul  ° 

(see  [8],  whore  the  procedure  for  the  solution  is  also  described) . 

The  time  interval  was;  selected  with  consideration  of  stability,  being 


r 
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equal  to: 


m.ix  : f v | i 


(1.18) 


The  three-point  formula  for  finding  the  line  perpendicular  to  the 
boundary  of  the  temperature  gradient  was  used  tc  calculate  the 
thermal  flux. 


The  calculations  were  irade  for  the  follcwirg  values  of  the 
parameters.  The  value  of  the  Franltl  number  was  fixed:  P = 5;  ‘1:p 
value  cf  the  Grashof  number  varied  up  to  G = 15«103;  the 
heterogeneity  parameter  was  assigned  two  values:  y = 0 and  y - 0.825. 


2.  Results  of  Calculations, 


The  main  results  are  shown  in  Figures  1-12  and  Tables  la  an  1 lb, 
where  the  integral  characteristics  of  convection  are  given  for  *he 
calculated  values  of  the  parameters.  All  of  the  integral 
characteristics  of  the  convective  process  were  considered  under 
stable  conditions.  Stable  ccnvect icn  conditions  were  achieved  ly  the 
transitional  process  of  attenuating  cscillaticrs  analogous  to  *he 
case  cf  homogeneous  viscosity  [9], 


Figures  2-3  give  the  dependences  (the  Erandtl  number  P = 5)  of 
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the  intensity  of  convective  movement  <|'<n  anc  the  dimensionless  heat 
flux  N cn  the  Grashcf  number  when  the  viscosity  is  constant  (y  = 0) 
and  there  is  a ten-fold  drcj  (y  - 1.637  for  the  square  and  r = 0.P25 
for  the  cylinder).  Figures  <a  and  3a  were  clotted  for  a cavity  w it  H a 
square  cross  section,  and  2t  and  3t  - for  a round  cylinder.  As  we  can 
see,  the  intensity  of  movement  and  heat  transfer  increase  somewha*  as 
a result  of  the  heterogeneity  cf  the  viscosity. 
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It  should  be  rioted  that  due  tc  the  asynmetry  of  the  t empe  re  tur‘~- 
field,  the  precision  of  the  computation  of  the  heat  fluxes  and 
turns  out  to  be  different  at  j i 0 (see  Tables  la  and  1b).  Figure  2 
shews  the  value  of  N depending  on  G,  equal  tc  ♦ N_)/2.  The 

convergence  of  dependences  "15(G)  at  large  Grsshcf  numbers  obviously 
indicates  an  increase  in  the  errer  cf  the  difference  system  comnared 
to  the  case  y = 0,  when  the  number  of  nonlinear  terms  in  the  equation 
cf  neve  me nt  is  much  smaller. 
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Table  1b. 


a 

1 -0 

500 

2.15 

1000 

2.68 

2500 

3.54 

5000 

! 4.21 

7500 

4.68 

10000 

1 5.02 

12  500 

5.31 

15000 

5.58 
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The  increase  in  the  intensity  of  mcvener.t  when  viscisi+v  ,'-per,s  j 
cn  temperature  is  illustrated  in  Fig.  4,  ir  which  the  maximum 
velocity  is  shown  in  dependence  on  the  htter  cgeneit  y parameter  y 
in  the  case  cf  a region  with  a sguare  cross  section.  In  the  range  of 
y in  question,  this  dependence  (G  = 2D»103,  P = 1 and  P = 5)  is  cl  osc 
linear. 


Fig.  4. 
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Me  will  point  out  that  the  integral  characteristics  of  the 
convective  process  vary  less  with  the  ircreese  in  y than  the  local 
characteristics.  In  the  case  of  a regicn  with  a sjuare  cross  section, 
ty  processing  the  results  fcr  p = c ami  G = 2C«103  we  obtain  tp-i 
following  af  proximate  dependences  ( 0 ,<  > ,<  1.6): 


^V(T)-/V(0)(l  +0.027),  (21) 

•V„(t)  = K(0)(1+0.02T). 

Figures  5-12  illustrate  *-hp  structure  cf  movement.  and  the 
temperature  field  at  different  values  of  the  parameters.  As  we  car 
see,  at  sufficiently  large  values  cf  tie  3ayleich  number  3a  = G*p, 
which  corresponds  to  Figures  5-f  and  Figures  11-12,  the  velocity  and 
temperature  boundary  layers  are  wel  1-develc  red  in  the  flow. 

Naturally,  the  effect  of  the  heterogeneity  cf  viscosity  increases 
with  tho  increase  in  paraseter  y.  This  effect  is  evident,  first,  in 
the  fact  that  flow  symmetry  increases  at.  y = C . The  flow  rates  on  t)e 
hot  wall  increase,  and  asymmetry  >f  fields  \p>  ard  T originates.  This 
asymmetry  is  evident,  ir  particular,  from  tr.e  shift  in  tne  isotherms, 
the  current  lines  and  tie  vertex  centers.  The  isotherm  corresponding 
tc  the  mean  temperature  (1  = 0)  shifts  toward  the  cold  fluid. 

Me  will  point  out  that,  the  patterns  of  movement  shown  in  Figures. 


9-12  qualitatively  agree  with  the  patterns  giver  in  experimental 
St- oeUj  L » o J • 
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Thus,  the  result'*  cf  the  calculations  indicate  that  the 
dependence  of  viscosity  on  temperature  causes  a more  or  less 
significant  change  in  the  local  ch a r acter is t ics  of  convective  flow. 
?s  for  the  integral  characteristics  - the  maximum  value  of  the 
current  function  and  the  dimensionless  heat  flux  - being  determined 
ly  the  mean  viscosity,  they  change  with  consideration  of  the 
temperature  dependence  cf  visccsi4-)  in  relatively  narrow  limits. 

The  authors  are  grateful  tc  G.  7.  Gershuni  and  Ye.  1. 
Zhukhcvitskiy  for  their  ccns+an4-  attention  tc  tie  work  and  their 
assistance  . 
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